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Results from unsteady Reynolds-averaged Navier-Stokes computations are described 
for two different synthetic jet flows issuing into a turbulent boundary layer crossflow 
through a circular orifice. In one case the jet effect is mostly contained within the bound- 
ary layer, while in the other case the jet effect extends beyond the boundary layer edge. 
Both cases have momentum flux ratios less than 2. Several numerical parameters are 
investigated, and some lessons learned regarding the CFD methods for computing these 
types of flow fields are summarized. Results in both cases are compared to experiment. 


Keywords: synthetic jet, Navier-Stokes, flow control, turbulence modeling 

I. Introduction 

Synthetic jets are attracting increasing interest in the aerospace community because they have 
many practical applications, including jet vectoring, control of separation, enhanced mixing, re- 
duction of wall skin friction, and virtual aeroshaping. As a practical tool, they are generally more 
attractive than steady blowing or suction because they do not require any complex fluid ducting or 
plumbing systems. Furthermore, synthetic jets have been found to achieve similar effectiveness to 
steady blowing or suction with considerably smaller momentum (and hence lower energy inputs). 1 

The problem of synthetic jet flow in a turbulent boundary layer has been studied both experi- 
mentally (see, for example 2-5 ) and computationally (see, for example 6-11 ). Earlier numerical work 
by Rizzetta et al. 12 for synthetic jets into quiescent flow demonstrated that accounting for the 
internal actuator geometry was important because it affected the jet profiles at the exit. Lee and 
Goldstein 6 used direct numerical simulation (DNS) to study an array of synthetic microjets in a 
crossflow. A 4-slot array was modeled, and adjacent jets were coupled 180° out of phase. As 
microjets, the Reynolds number of the jets was very small, on the order of 10. Weak jets were 
confined to the linear sublayer, but when the jets penetrated up to the buffer region, they affected 
the events of the boundary layer significantly. 

Mittal et al. 7 and Mittal and Rampunggoon 8 used 2-D Navier-Stokes to investigate the use 
of stronger synthetic jets immersed in a boundary layer to form large mean recirculation bubbles 
(for aeroshaping). They modeled the cavity, but their studies were purely numerical (there was no 
associated experiment). Similarly, Ravi et al. 9 performed a purely numerical study of synthetic jet 
in crossflow. They also modeled the cavity, and used 3-D Navier-Stokes to investigate the effects 
of a square and two rectangular shaped slots. 

Cui et al. 10 performed 2-D computations on a synthetic jet in a turbulent crossflow bound- 
ary layer using unsteady Reynolds-averaged Navier-Stokes (URANS) with the Spalart-Allmaras 
turbulence model. 13 They compared their computational fluid dynamics (CFD) results with two 
experiments: one in which the discrete vortices scaled with the boundary layer thickness, and 
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one in which they scaled with the inner viscous sublayer. Overall, they found qualitatively good 
agreement with experiment, although in the latter case the computed vortices moved too slowly 
downstream. Also, in both cases there were significant quantitative differences between the CFD 
results and experiment. 

A CFD validation workshop on synthetic jets and turbulent separation control (CFDVAL2004) 14 - 15 
was held in March 2004 in Williamsburg, VA. The purpose of the workshop was to assess the cur- 
rent capabilities of different classes of turbulent flow solution methodologies to predict flow fields 
induced by synthetic jets and separation control. Case 2 of the workshop was a synthetic jet is- 
suing into a turbulent boundary layer through a circular orifice in the floor. This experiment was 
conducted at NASA Langley Research Center specifically for use as a CFD validation case. 16 A 
similar experiment was also recently conducted at NASA Glenn Research Center. 17 

A summary of the CFD results for case 2 of the CFDVAL2004 workshop is given in Rumsey 
et al. 15 There were five contributers who ran ten separate cases. One contributor used large eddy 
simulation (LES) and the others all used URANS. On the whole, reasonably good qualitative re- 
sults were obtained compared with experiment both for long-time averaged and phase averaged 
quantities, but there were significant variations among the CFD results. Different turbulence mod- 
els were found to have a much smaller effect than different grids, codes, and boundary conditions. 
Cui and Agarwal 11 also computed case 2 from CFDVAL2004, using both URANS with Menter’s 
SST turbulence model 18 and detached eddy simulation (DES). 19 The two methods produced over- 
all similar results in qualitative agreement with experiment, and DES was able to simulate some of 
the larger cross-stream velocity components downstream. 

The intent of the current paper is to summarize, compare, and contrast the CFD results obtained 
using the URANS code CFL3D 20 for the two cases hereafter described as the “CFDVAL case” 
(Schaeffler and Jenkins 16 ) and the “NASA Glenn case” (Milanovic et al. 17 ). Previous CFD results 
using CFL3D were reported for the CFDVAL workshop case, 21,22 but only limited results were 
reported for the NASA Glenn case in a paper devoted primarily to describing the experiment. 17 
Furthermore, a summary is given of lessons learned regarding the CFD methods for computing 
these types of flow fields, including grid size, time step, subiterations, boundary conditions, and 
turbulence modeling. 


II. Numerical Method 

The computer code CFL3D 20 solves the three-dimensional, time-dependent, Reynolds-averaged 
compressible Navier-Stokes equations with an upwind finite-volume formulation (it can also be ex- 
ercised in two-dimensional mode of operation for 2-D cases). Because the code is compressible, 
the Navier-Stokes equations are averaged using Favre averaged variables; e.g., fm,/ p. CFL3D 
can solve flows over multiple-zone grids that are connected in a one-to-one, patched, or overset 
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manner, and can employ grid sequencing, multigrid, and local time stepping when accelerating 
convergence to steady state. Upwind-biased spatial differencing is used for the inviscid terms, and 
flux limiting is used to obtain smooth solutions in the vicinity of shock waves, when present. Vis- 
cous terms are centrally differenced, and cross-diffusion terms are neglected. For very low Mach 
number flows, preconditioning 23 is used to insure convergence and accuracy of the solutions. 

The CFL3D code is advanced in time with an implicit approximate factorization method. The 
implicit derivatives are written as spatially first-order accurate, which results in block tridiagonal 
inversions for each sweep. However, for solutions that utilize Roe flux-difference splitting, 24 the 
block tridiagonal inversions are further simplified using a diagonal algorithm with a spectral radius 
scaling of the viscous terms. 

In time-accurate mode, CFL3D uses pseudo-time stepping with multigrid and achieves second 
order temporal accuracy. With pseudo-time stepping, subiterations are used to reduce the lineariza- 
tion and factorization errors, and advance the solution in pseudo-time to the next physical time. For 
a non-deforming mesh, the time dependent compressible Navier-Stokes equations can be written 
as 


J dt 


R{ Q) 


(i) 


where Q is the vector of conserved variables \p, pu, pv, pw, e] T , J is the Jacobian of the generalized 
coordinate transformation, and R is the right-hand- side: 
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made up of gradients of the inviscid and viscous flux terms. Eq. (1) can be discretized with back- 
ward differencing: 
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The superscripts n, n + 1, and n — 1 indicate the time level. When 0 = 0 the method is first- 
order temporally accurate, and when 0 = 1/2 the method is second-order temporally accurate. To 
employ pseudo-time stepping, a pseudo time term is added to Eq. (3), and the equation is iterated 
in m, where m is the subiteration counter. After linearizing R and performing some additional 
manipulation, one obtains: 
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where r is the pseudo time parameter, AQ m = Q m+1 — Q m , and 
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Eq. (4) is approximately factored and written in primitive variable form; it is solved as a series of 
sweeps in each coordinate direction. Additional details are given in Krist et al. 20 In the current 
study, the effects of different physical time steps and number of subiterations were investigated. 

The turbulence models are solved uncoupled from the mean flow equations using implicit ap- 
proximate factorization. Their advective terms are solved using first-order upwind differencing. 
Many turbulence models are available in CFL3D, but only those used in the current study are 
mentioned here. Descriptions of the one-equation Spalart-Allmaras (SA) and the two-equation 
Menter k-u shear-stress transport (SST) turbulence models can be found in their respective refer- 
ences. 13,18 Both of these models are linear eddy-viscosity models that make use of the Boussinesq 
eddy-viscosity hypothesis. The EASM model is a nonlinear eddy-viscosity model in k-u form: 
it stands for explicit algebraic stress model. A detailed description of this model can be found in 
Rumsey and Gatski. 25 

For time-accurate URANS computations, it can be shown that the dependent variables calcu- 
lated for forced periodic cases correspond to phase- averaged variables. Following Reynolds and 
Hussain, 26 triple decomposition can be used to represent a fluctuating signal: 


/ ( x > t) = f (x) + / (x, t) + /' (x, t) 


(6) 


where / (x) is the (time) mean value, / (x, t) is the statistical contribution of the organized motion, 
and f (x, t) is the turbulence. The long-time average is defined as 

/( x )= lim 4/ f{x,t)dt (7) 

T^oo 1 Jo 

and the phase average is defined as 


1 N 

(/ ( x ? t)) = lim — / (x, t + nr) (8) 

N ^°° N n = o 

where r in the case of a synthetic jet flow field is the period of the imposed oscillatory motion. 
Thus, the wave component is 

/Cm) = (/( x A)) - /( x ) (9) 

and the decomposition of Eq. 6 can be written as 


/( x ,t) = (/ ( X A)) + /' ( x ,t)- 


( 10 ) 
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This form now resembles standard Reynolds splitting, except that the flow is split into a slowly 
varying mean (phase average) and a random fluctuating part. See also Gatski and Liu. 27 As a result, 
the final conservation equations in terms of phase-averaged variables ( ) are identical in form to the 
standard Reynolds-averaged equations. In other words, when URANS is used for a forced periodic 
flow field such as that produced by a synthetic jet, the computation typically eventually settles 
down to a nearly exactly repeatable periodic variation. If and when this repeatability is attained, 
any single point during the cycle corresponds to a phase-averaged result from the experiment. An 
example showing the cycle-to-cycle repeatability for a global integrated quantity is shown in Fig. 1 
for the CFDVAL case. Here, coefficient of drag nondimensionalized by a unit reference length 
(the absolute level of which is not particularly meaningful in this case) is plotted as a function of 
iteration. An FFT analysis (not shown) over the last three periods confirms that repeatability has 
been obtained, with the coefficients of lower frequency modes more than 2 orders of magnitude 
lower than that of the driving frequency. 

All of the results for the current two cases were expected to be symmetric about the center 
plane of the orifice. The CFDVAL case was computed with a full-plane grid, and results indeed 
maintained symmetry as expected. For the NASA Glenn case, use of a full-plane grid caused 
asymmetries to develop while establishing an initial steady-state flow field with constant blowing. 
The cause of these asymmetries is not precisely known, but it appears to be related to the strength of 
the blowing. The jet in the NASA Glenn case was strong enough to penetrate through the boundary 
layer to the freestream. If the jet strength was lowered significantly, then the full-plane constant 
blowing flow field remained symmetric. In any case, even with asymmetric initial conditions, 
these asymmetries went away and the solution became symmetric after the cyclic time-dependent 
boundary condition, to be described in the next section, was initiated. 

III. Results 

For the CFDVAL case, flow passed in and out of a circular orifice D = 6.35 mm in diameter. 
The orifice was located on the floor of a wind tunnel splitter plate with a turbulent boundary layer 
atM=0.1 (Uoo = 34.6 m/s) and approximate boundary layer thickness of21 mm ( 5/D = 3.3). The 
Reynolds number was 2230 per mm, or 14,160 per orifice diameter. The jet was driven electro- 
mechanically by a bottom-mounted square- shaped rigid piston mounted on an elastic membrane 
inside the cavity chamber beneath the splitter plate. The cavity was approximately 1.7 mm deep 
on average while the piston was operating, and the piston moved approximately ±0.77 mm. The 
thin volume inside the cavity varied during the operation of the piston, from roughly 1.0 x 10 4 - 
2.6 x 10 4 mm 3 . The frequency was 150 Hz, and the maximum velocity (w max ) out of the orifice 
was approximately 1. 311^. The momentum flux ratio, defined as (u' m ».x/U^) 2 , was approximately 
1.69. In this case, the effects due to the jet were mostly contained within the boundary layer. 
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Experimental data were measured using laser-doppler velocimetry (LDV) and particle image ve- 
locimetry (PIV). 16 

For the NASA Glenn case, flow passed in and out of a circular orifice D = 19.05 mm in 
diameter. The Mach number of the freestream was M=0.0175 (U^ = 6.1 m/s), and the turbulent 
boundary layer thickness was approximately 7.6 mm (5/D = 0.4). The Reynolds number was 
395.3 per mm, or 7530 per orifice diameter. The jet was driven by a loud speaker at the bottom 
of a cavity mounted beneath the tunnel floor. The cavity was 152 mm deep (with a volume of 
approximately 2.8 x 10 7 mm 3 ). The frequency was 24 Hz, the maximum velocity out of the orifice 
was approximately l.Of/oo, and the momentum flux ratio was approximately 1.0. In this case the 
effects due to the jet extended well beyond the edge of the boundary layer. Experimental data were 
measured using hot-wire anemometry. 17 

For all of the results to follow, the ^-direction is downstream and ^-direction is up. The y = 0 
plane represents the center plane aligned with the center of the orifice. Note that for the two 
cases the phases (which were designated by the individual experiments) were not the same: for 
the CFDVAL case peak discharge occurred near phase = 120° and peak suction near phase = 300°, 
whereas for the NASA Glenn case peak discharge occurred near phase = 260° and peak suction 
near phase = 80°. However, this phase designation is arbitrary. 

Two different gridding approaches were taken for modeling the two cases. Because the CFD- 
VAL case employed such a shallow cavity, the computational grid included the (as-built) cavity. 
A time-dependent velocity boundary condition was applied at the (fixed) bottom of the cavity to 
simulate the effect of a moving wall. Additional details can be found in Rumsey. 22 For the NASA 
Glenn case, the cavity was so large relative to the orifice that it was impractical to model it in the 
computations. Therefore, in this case the grid included the orifice only, which extended 25.4 mm 
below the tunnel floor. For both cases, the initial conditions were steady-state flow with constant 
blowing. Subsequently, a time-dependent velocity boundary condition was applied at the bottom 
of the orifice or cavity. This boundary condition was calibrated in order to attempt to match the 
velocity as a function of phase at the orifice exit, as measured experimentally. The time-dependent 
velocities prescribed were: 

U — 0 V = 0 W = [ (pW ) max / p\ COS ( 2 TfFt) (11) 

where (pw) max = O.OOOSpooCioo for the CFDVAL case, and 0.0175poo<7oo for the NASA Glenn case. 
F was the frequency of oscillation. Due to the very low freestream Mach number of the NASA 
Glenn case, preconditioning was employed in that particular computation. At the upstream grid 
boundary, a turbulent profile was specified that matched the corresponding experimental incoming 
boundary layer profile. 

Figure 2 shows time history of the velocities near the center of the orifice for the CFDVAL 
case. Results are shown for two different full-plane grids (the fine grid had 4. 1 million total points 
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and the medium grid was created using every other point from the fine grid), and for three dif- 
ferent turbulence models (SA=Spalart-Allmaras, SST=Menter’s shear stress transport k-u, and 
EASM=explicit algebraic stress model). For the most part differences due to grid and turbulence 
model were relatively small. The CFD mimicked the general overall trends measured in the exper- 
iment, but there were some noticeable specific differences. However, these results were the best 
it was possible to attain within the constraints imposed by the modeling and boundary conditions 
used. 

Figure 3 shows long-time average //-velocity contours for the CFDVAF case on the center plane 
using CFD with the SA model. Although not shown, results from the other turbulence models 
looked very similar. The mean perturbation to the velocity field was mostly contained within the 
boundary layer. There were not enough experimental data available to cover the same plotting 
range as Fig. 3, but Fig. 4 shows a comparison between CFD and experiment (PIV) for //-velocity 
contours in a zoomed-in region very near the orifice exit, where xq is the orifice center. There 
was only a crude qualitative correspondence between the two. The CFD result predicted a larger 
mean perturbation to the flow field (including a very small reverse flow region near the orifice lip at 
(x — xo )/D = 0.5 that is difficult to see in the figure), whereas the experimental PIV data showed 
the retarded velocity region to be smaller and offset downstream from the lip. 

CFD results compared better with the FDV experimental data for the CFDVAF case. A com- 
parison of w-velocity profiles at phase = 120° (at (x — x 0 )/D = 1.0) is shown in Fig. 5. All three 
turbulence models yielded similar results in good agreement with the data at this critical phase 
when the boundary layer was being severely perturbed by the vortex created during the expulsion 
part of the synthetic jet cycle. The figure also shows a comparison between the FDV and PIV data: 
there were fairly significant differences, the cause of which is unknown. Overall for the CFD- 
VAF case, CFD was generally successful in qualitatively predicting many flow features - both 
long-time-averaged and phase-averaged - compared with FDV measurements. 

For the NASA Glenn case, only one turbulence model, SA, was employed. As discussed in 
section II, there were some problems maintaining symmetry on a full-plane grid when establishing 
the constant blowing initial condition, but subsequent unsteady computations re-established sym- 
metry. The effect of computing with a full-plane grid vs. a half-plane grid was investigated for this 
case. A sample comparison showing long-time average //-velocity contours in a plane parallel to 
the floor is shown in Fig. 6. Results showed very little difference, so all NASA Glenn results for 
the remainder of the paper were obtained using the half-plane grid only. 

Figure 7 shows the fine half-plane grid, made up of 1.8 million points. This grid had the same 
(rescaled) distribution of grid points as the CFDVAF case outside of the orifice. A medium grid 
was also employed that was created using every other point from the fine grid. The time history of 
velocity at the orifice center is shown in Fig. 8. Two different grid sizes made very little difference 
in the results. The total velocity agreed fairly well with the experiment, which did not measure 
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the individual velocity components like in the CFDVAL case. The largest discrepancy between 
CFD and experiment was in the magnitude during the suction part of the cycle. The computations 
accurately captured the double hump in total velocity during the discharge part of the cycle seen in 
the experiment. 

These results were obtained using a time step corresponding to 720 steps per cycle, and 10 
subiterations per time step. The effects of time step and number of subiterations were also investi- 
gated for the NASA Glenn case. Total velocity at the orifice center is plotted in Fig. 9 for the three 
conditions of 720 steps per cycle and 10 subiterations, 720 steps per cycle and 20 subiterations, 
and 360 steps per cycle and 10 subiterations. These results were obtained on the medium grid. 
There were imperceptible differences between the three results. Contour plots in the flow field 
also indicated a high degree of temporal convergence. For example, Fig. 10 shows phase-averaged 
u - velocity contours in the plane z = 0.192 mm at phase = 0°. The left half of the figure shows 
results using 360 steps per period and the right half shows results using 720 steps per period. 

Figure 1 1 shows long-time average contours of u/U , ^ in the center plane for the NASA Glenn 
case. The agreement between CFD and experiment was qualitatively very good. The mean per- 
turbation to the velocity field in this case extended well beyond the boundary layer edge. Phase 
averaged contour plots also showed good qualitative agreement. For example, Fig. 12 compares 
streamwise velocity contours from CFD and experiment at phase = 240° in the center plane. In 
particular, CFD captured the size and position of the disturbance very well as it convected down 
the plate after being emitted from the orifice, although it predicted lower velocity in the low- 
momentum fluid and the “tilt” in the structure predicted by the CFD was less pronounced than in 
the experiment. 

Quantitative comparisons of w-velocity profiles as a function of phase at a fixed position (x — 
x 0 )/D = 2.235, y — 0 are shown in Fig. 13. From Milanovic et al., 17 the convection speed of 
the synthetic jet structure created during the peak expulsion part of the cycle (phase ~ 260°) was 
approximately 0.75f/oo. Thus, the structure convected to (x — x 0 )/D = 2.235 at approximately 
phase = 340°. As the figure shows, at the two phases 320° and 0° surrounding this “event,” the 
velocity profiles were significantly affected, and the CFD predicted perturbed velocities in the 
boundary layer that were too low compared with experiment. However, the profile shapes were 
qualitatively similar, and the profiles at the other three phases were quantitatively predicted very 
well. The long-time average also yielded excellent agreement. Figure 13 also compares results 
on both the fine and medium grids, giving an indication of the approximate effect of discretization 
error for these grid sizes. The effect was relatively small. 

Figures 14 and 15 are contour plots showing the vertical velocity at the orifice for both the 
CFDVAL and NASA Glenn cases, within 20° of peak discharge and peak suction. In the figure 
the closest phase available from the stored CFD data is displayed (recall that the arbitrary phases 
between the two experiments were different). Both cases showed similar behavior. Near peak 
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discharge, the velocities out of the orifice were very nonuniform, and the strongest velocities oc- 
curred near the side edges of the orifice, downstream of center. Near peak suction, the velocities 
were more uniform over much of the orifice area. These figures indicate the importance of mod- 
eling the inside of the orifice in CFD computations. Although a simple time-dependent top-hat 
profile might adequately model the relatively uniform inflow during the suction phase, it would 
be difficult to prescribe a time-dependent profile at the exit plane that mimicked the non-uniform 
velocities present during the expulsion part of the cycle. The effect of ignoring the inside of the 
orifice, and instead applying the time dependent boundary condition (Eq. 1 1) directly at the orifice 
exit, was investigated. Quantitative comparisons are given in Fig. 16, showing both phase-average 
results as well as the long-time average. Not modeling the inside of the orifice had an overall 
negative impact on the comparison with experiment, as expected. 

IV. Conclusions 

Time-dependent Reynolds-averaged Navier-Stokes computations were compared for two dif- 
ferent synthetic jet flows issuing into a turbulent boundary layer crossflow through a circular ori- 
fice. In one case the jet effect was mostly contained within the boundary layer, while in the other 
case the jet effect extended beyond the boundary layer edge. The former case was originally one 
of the subjects of the CFDVAL2004 workshop on synthetic jets. 

The use of URANS was shown to be appropriate and viable for synthetic jet flows. In both 
cases described here, the computations did a reasonably good job predicting the qualitative behav- 
ior of the synthetic jet flow field. However, quantitatively there were many specific differences. 
Several numerical parameters were investigated, including effect of grid size, time step, number 
of subiterations, and turbulence model. Differences between results on the fine and medium grids 
were relatively small. Using at least 360 time steps per period and at least 10 subiterations per 
time step was more than sufficient to produce results that did not change perceptibly with further 
temporal refinement. Different turbulence models for the CFDVAF case had some influence on the 
solution, but these differences were generally not too large. 

For synthetic jets like those considered here, either a full-plane grid or a half-plane grid (with 
symmetry) could be used. The current computations showed the importance of including the region 
inside the orifice when modeling the problem. Use of time-dependent top-hat boundary conditions 
at the orifice exit plane was an oversimplification that failed to capture the complex nature of the 
flow field near the orifice, particularly during the expulsion phase of the cycle. 
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Figure Legends: 


Fig 1. Scaled drag coefficient as a function of iteration, showing periodic repeatability of 
URANS computation for CFDVAL case. 

Fig 2. Time history of velocity near center of orifice for CFDVAL case (z = 0.4 mm above the 
floor). 

Fig 3. Computed long-time-average u/U^ contours in center plane of CFDVAL case (orifice 
is centered at (x — x 0 )/D = 0), SA, fine grid. 

Fig 4. Close-up of long-time-average u/U^ contours in center plane of CFDVAL case near 
orifice. 

Fig 5. Profiles of w-velocity at (x — x 0 )/D = 1.0, phase = 120°. 

Fig 6. Contours of long-time average ^-velocity in the z = 0.192 mm plane for NASA Glenn 
case using SA, medium grids, showing difference between full-plane computations (left half) and 
half-plane with symmetry (right half). 

Fig 7. Picture of the half-plane 6-zone grid used for NASA Glenn case (1.8 million grid points 
total). 

Fig 8. Time history of velocity at center of orifice for NASA Glenn case, SA. 

Fig 9. Effect of time step and number of subiterations on time history of total velocity at center 
of orifice for NASA Glenn case, SA, medium grid. 

Fig 10. Contours of ^-velocity at phase = 0° in the z = 0.192 mm plane for NASA Glenn case 
using SA, medium grid, showing difference between 360 steps per period (left half) and 720 steps 
per period (right half). 

Fig 11. Comparison of long-time average u/U^ contours in center plane of NASA Glenn case 
(orifice is centered at (x — x 0 )/D = 0), SA, fine grid. 

Fig 12. Comparison of phase-averaged u/U^ contours in center plane of NASA Glenn case 
(orifice is centered at (x — x 0 )/D = 0), phase = 240°, SA, fine grid. 

Fig 13. Profiles of w-velocity at (x — x 0 )/D = 2.235, y — 0 as a function of phase and 
long-time average for NASA Glenn case, SA. 

Fig 14. Contours of u’-vclocity at orifice exit plane near peak discharge part of cycle using SA, 
fine grid; CFDVAL case on left, NASA Glenn case on right. 

Fig 15. Contours of w-velocity at orifice exit plane near peak suction part of cycle using SA, 
fine grid; CFDVAL case on left, NASA Glenn case on right. 

Fig 16. Effect of orifice treatment on phase-average and long-time-average //-velocity profiles 
at (x — x q )/D = 2.235, y = 0 for NASA Glenn case, SA, medium grid. 
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Figure 2. Time history of velocity near center of orifice for CFDVAL case (z = 0.4 mm above the floor). 
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(a) comparison of CFD with experiment (LDV) (b) comparison of LDV with PIV measurements 




Figure 5. Profiles of u-velocity at (x — xo) / D = 1.0, phase = 120°. 
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Figure 6. Contours of long-time average (/-velocity in the z = 0.192 mm plane for NASA Glenn case using SA, 
medium grids, showing difference between full-plane computations (left half) and half-plane with symmetry 
(right half). 
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Figure 7. Picture of the half-plane 6-zone grid used for NASA Glenn case (1.8 million grid points total). 
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Figure 8. Time history of velocity at center of orifice for NASA Glenn case, SA. 
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Figure 9. Effect of time step and number of subiterations on time history of total velocity at center of orifice 
for NASA Glenn case, SA, medium grid. 
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Figure 10. Contours of (/-velocity at phase = 0° in the z = 0.192 mm plane for NASA Glenn case using SA, 
medium grid, showing difference between 360 steps per period (left half) and 720 steps per period (right half). 
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Figure 11. Comparison of long-time average u/U^ contours in center plane of NASA Glenn case (orifice is 
centered at ( x - x 0 )/D = 0), SA, fine grid. 


23 of 27 



( x *Xq)/D 



( x *Xq)/D 


Figure 12. Comparison of phase-averaged u/Uoo contours in center plane of NASA Glenn case (orifice is cen- 
tered at (x - x 0 )/D = 0), phase = 240°, SA, fine grid. 
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(a) CFDVAL case 

w/Uinf contours at phase=1 20, z=0 



(b) NASA Glenn case 

w/Uinf contours at phase=240, z=0 



Figure 14. Contours of w-velocitv at orifice exit plane near peak discharge part of cycle using SA, fine grid; 
CFDVAL case on left, NASA Glenn case on right. 


(a) CFDVAL case 

w/Uinf contours at phase=280, z=0 



(b) NASA Glenn case 

w/Uinf contours at phase=80, z=0 



Figure 15. Contours of w-velocity at orifice exit plane near peak suction part of cycle using SA, fine grid; 
CFDVAL case on left, NASA Glenn case on right. 
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Figure 16. Effect of orifice treatment on phase-average and long-time-average (/-velocity profiles at ( x 
x 0 )/D = 2.235, y = 0 for NASA Glenn case, SA, medium grid. 
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